Hysteresis effects in the phase diagram of multiferroic GdMnOa 
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We present high-resolution thermal expansion a(T) and magnetostriction AL(H)/L measure- 
ments of GdMnC>3, which develops an incommensurate antiferromagnetic order (ICAFM) below 
Tn — 42 K and transforms into a canted A-type antiferromagnet (cAFM) below T c ~ 20 K. In addi- 
tion, a ferroelectric polarization P||a is observed below Tfe for finite magnetic fields applied along 
the b direction. In zero magnetic field we find a strongly anisotropic thermal expansion with certain, 
rather broad anomalous features. In finite magnetic fields, however, very strong anomalies arise at 
T c for fields applied along each of the orthorhombic axes and at Tfe for fields along the b axis. 
Both phase transitions are of first-order type and strongly hysteretic. We observe a down-bending 
of the ICAFM-to-cAFM phase boundary T C (H) for low magnetic fields and our data give evidence 
for coexisting phases in the low-field low-temperature range. 



The recent discovery of very large magnetoelectric ef- 
fects in the rare-earth manganites _RMn03 has reopened 
the field of the so-called multiferroic materialsi^ Mul- 
tiferroic means, that several ferro-type orders like fer- 
romagnetism, ferroelectricity or ferroelastivity coexist. 
The rare-earth manganites -RMnC>3 may be grouped into 
the hexagonal ones (R— Ho,. . . , Lu) and those with 
orthorhombically distorted perovskite structures (R= 
La,. . . , Dy). The hexagonal RMn03 show both, ferro- 
electric and magnetic order, but the respective ordering 
temperatures differ by an order of magnitude, fpE — 
1000 K and T N ~ 100 In contrast, RMn0 3 with R = 
Gd, Tb, and Dy possess comparable transition temper- 
atures for the magnetic and the ferroelectric orderingi^ 
Both, TbMnC>3 and DyMn03 develop an incommensu- 
rate antiferromagnetic order (ICAFM) below about 40 K. 
The incommensurability continuously changes upon cool- 
ing and becomes almost constant below T c ~ 28 K and 
~ 19 K for R = Tb and Dy, respectively. The transition 
to this long-wavelength incommensurate antiferromag- 
netic (LT-ICAFM) phase is accompanied by ferroelectric 
(FE) ordering with a polarization P||c, which flops to 
P| | a above a critical magnetic field applied along the a or 
b direction, while it is suppressed for large fields along a 4 
GdMnO-3 also shows an ICAFM order below T N ~ 42 K 
and a second transition at T c ~ 23 K. Based on the 
observed weak ferromagnetism^ and on X-ray diffrac- 
tion studies^ a canted A-type antiferromagnetic ordering 
(cAFM) has been proposed for T < T c , but a direct mag- 
netic structure determination has not yet been published. 
Both transition temperatures weakly increase in a mag- 
netic fields Concerning the FE polarization, contradic- 
tory results have been reported. Kuwahara et alm> find a 
finite polarization below 13 K, while Kimura et alX ob- 
serve FE order only between 5 K and 8 K. The magnitude 
of P is much smaller than in TbMnOs and DyMnC>3 and 
the direction is P||a. Moreover, there are no magnetic- 
field-induced polarization flops. Instead the polarization 
is stabilized by a magnetic field along 6, while it is im- 
mediately suppressed for fields along a and 



In order to study the coupling of the various phase 
boundaries to lattice degrees of freedom we have con- 
ducted high-resolution measurements of thermal expan- 
sion and magnetostriction on GdMnOs. We find pro- 
nounced anomalies at all transitions, i.e. all transitions 
strongly couple to the lattice. Thus, our data allow for 
a precise determination of the magnetic-field tempera- 
ture phase diagram and clearly reveal that the ICAFM- 
to-cAFM as well as the FE transition are of first oder 
with strong hysteresis. For low magnetic fields we find 
a down-bending of the ICAFM-to-cAFM phase bound- 
ary and evidence for a coexistence of both phases in this 
low-field range. This should be taken into account for 
further experiments, in particular for zero magnetic field. 
Probably, this may also explain the contradictory results 
reported for GdMnOai&S 

The GdMnOa single crystal used in this study is a 
cuboid of dimensions 1.7x2x1. 45 mm 3 along the a, b, and 
c direction (Pbnm setting), respectively. It was cut from 
a larger crystal grown by floating-zone melting. Phase 
purity was checked by x-ray powder diffractionsiS Mag- 
netization, resistivity and specific heat data of the same 
crystal are reported in Ref. @- The linear thermal expan- 
sion aj = dlnLi/ dT and magnetostriction ALi (if) /Li = 
[Li (if) — Lj(0)]/Lj(0) have been measured by a home- 
built high-resolution capacitance dilatometer^i Here, Li 
denote the lengths parallel to the different crystal axes 
i = a, b, and c. In general, we studied the length changes 
in longitudinal magnetic fields up to 14 T, i.e. H\\i, and, 
in addition, we measured a.\> for H\\c. 

Fig- HI shows the zero-field a, of GdMnOs for i = a, b, 
and c, which are strongly anisotropic and show several 
anomalies. The sharp anomalies around 41 K signal the 
Neel transition of the Mn ions and their shape is typ- 
ical for a second-order phase transition. According to 
previous publications further anomalies are expected at 
lower T: (i) at the ICAFM-to-cAFM transition around 
T c ~ 23 K, (ii) around Tfe ^ 10 K, where the FE order- 
ing sets in, and (iii) around 6 K due to the ordering of 
the Gd moments. Indeed, there are pronounced anoma- 
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FIG. 1: (Color online) Thermal expansion of GdMn03 along 
a, b, and c. 

lies around 6 K, small ones around 10 K, and interme- 
diate ones around 20 K with different signs and magni- 
tudes for the different i. However, all these anomalies 
are rather broad making a clear identification of transi- 
tion temperatures difficult. This drastically changes for 
finite magnetic fields. 

In Fig. [5] we show a a (T) for H\\a. With increasing 
magnetic field a broad anomaly shows up around 10 K, 
changes sign and smears out above 2 T. The origin of 
this anomaly is unclear, it may be related to the com- 
plex interplay of the magnetism of the Mn and Gd ions4i 
The most drastic change occurs, however, around 18 K, 
where a huge anomaly emerges between and 1 T. This 
anomaly is close to the observed ICAFM-to-cAFM tran- 
sition at T c and systematically shifts to higher T with 
further increasing field, in agreement with the observed 
field dependence T c (i?)iii The negative sign of the 
anomaly means that the transition from the cAFM to 
the ICAFM phase, with increasing T, is accompanied by 
a pronounced contraction of the a axis (see also Fig. [2J . 
The a a (T) curves of Fig. [2^, have been recorded with in- 
creasing T in the so-called field-cooled (FC) mode, i.e. 
the field has been applied at T ~ 50 K. We have also 
recorded the data during the cooling runs. Fig. ^jp com- 
pares a a {T) obtained with increasing and decreasing T 
for H = 4 T. Obviously, there is a strong hysteresis at the 
ICAFM-to-cAFM transition identifying this transition as 
first-order type. 

The lower panels of Fig.[3]present the magnetostriction 
AL a (H)/L a for H\\a at constant T. At 20 K an anoma- 
lous expansion of the a axis occurs around 4 T, which 
is related to a transition from the ICAFM to the cAFM 
phase as a function of increasing H. Both, the position 
of the AL a (H)/ L a anomaly in the H — T plane and its 
magnitude fit to the position and magnitude of the cor- 
responding a a {T) anomaly due to the ICAFM-to-cAFM 
transition as a function of decreasing T. The AL a (H) jL a 
curves obtained with increasing and decreasing H also 
show a hysteresis at the ICAFM-to-cAFM transition. 
With decreasing T the anomaly of AL a (H)/ L a shifts to- 
wards lower field. Around 12 K an anomalous expan- 
sion of the a axis is observed with increasing H, how- 
ever, this expansion is not reversed upon decreasing the 
field. Due to the large hysteresis GdMn03 remains in 
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FIG. 2: (Color online) Panel (a): Field-cooled thermal expan- 
sion a a (T) for H\\a measured with increasing T. Panel (b) 
shows the hysteresis around T c between the data obtained 
with increasing and decreasing T. Panel (c) compares two 
zero-field a a (T) curves for increasing T. The first curve [0T- 
(1)] has been recorded after zero-field cooling the sample and 
the second [0T-(2)] after cooling the sample to 12 K and sub- 
sequently applying and removing a magnetic field of 8 T. The 
lower panels (d) and (e) present the zero-field cooled magne- 
tostriction ALi(H)/Li recorded with increasing (solid lines) 
and decreasing (dashed) field. 



the cAFM phase after the magnetic field is switched off. 
With further decreasing T, the anomaly of AL a (H) / L a 
again shifts to higher field and at 5 K the anomalous 
expansion occurring with increasing H is reversed again 
with decreasing field. Apparently, we have traced the 
field-induced transition to the cAFM phase down to our 
lowest T. This means that the ICAFM-to-cAFM phase 
boundary T C (H) shows a clear down-bending in the low- 
field range (see Fig. [5J. In order to verify this, we have 
carried out the following zero-field a a (T) measurement: 
after cooling the sample to 12 K in zero field we have ap- 
plied a magnetic field of 8 T in order to enter the cAFM 
phase and due to the hysteresis of the ICAFM-to-cAFM 
transition the sample should remain in the cAFM phase 
after removing the field again. Fig. [2t shows a compar- 
ison of this zero-field a a (T) [labeled as 0T-(2)] with the 
conventional FC a a (T) [0T-(1)]. In contrast to the 0T- 
(1) curve, the 0T-(2) measurement shows a sharp peak 
around T c , which is comparable to the anomaly of the 
finite-field ot a (T) curves. This clearly confirms the down 
bending of the ICAFM-to-cAFM phase boundary. Thus, 
the (pure) cAFM phase of GdMn03 cannot be reached 
by cooling the crystal in zero field. The broad anomalies 
in the zero-field ati(T) curves (see Fig.^l probably arise 
from a partial ICAFM-to-cAFM transition and suggest 
that both phases coexist in the low-field region. 
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FIG. 3: (Color online) The relative length changes of the a, b, 
and c axes as a function of T at constant field (H = 6 T; left) 
and as a function of H at constant temperature (T = 15 K; 
right). In all cases H was applied parallel to the measured 
direction. 



In Fig. we compare the relative length changes as a 
function of field and temperature for all three axes. The 
c axis behaves similar to the a axis. The anomalies of 
c have the same signs but smaller magnitudes as com- 
pared to those of a and, moreover, their position for the 
same field (temperature) is located at higher T (lower 
H). Since in all cases we applied H\\Li, the latter differ- 
ence signals the anisotropy with respect to the field direc- 
tion. Our results (not shown) for a c (T) and AL C (H)/ L c 
for other fields and temperatures, respectively, are very 
similar to those obtained for the a axis. In particular, 
the distinct expansion at the ICAFM-to-cAFM transi- 
tion as a function of field or temperature is only present 
for H > 0.2 T. The hysteresis is also strong for H\\c, but 
less pronounced than for H\\a (see Fig. |SJ). A different 
phenomenology is observed for ALb/Lb for H\\b. Firstly, 
the anomalies at the ICAFM-to-cAFM transition are of 
opposite signs and, secondly, an additional anomaly oc- 
curs at a lower temperature. 

In Fig. 0| we present otb(T) for different H\\b. Again, 
a very pronounced anomaly evolves around T c ~ 18 K 
in small fields and shifts to higher T with further in- 
creasing field. The behavior of a& is analogous to our 
results on a a and a c , only the signs of the anomalies are 
different. In contrast to a a and a c , however, an addi- 
tional, very pronounced anomaly develops around 12 K 
in finite fields. Based on the polarization data of Ref. 0, 
we attribute this anomaly to the FE ordering at Tfe- 
As a further verification we have also measured ab in a 
transverse magnetic field H\\c and did not find such an 
additional anomaly (see inset of Fig. 0Ji). As shown in 
the lower inset of Fig. 0Jd the transition to the FE phase 
is also of first-order type with a broad hysteresis. With 
further increase of the magnetic field above 8 T the mag- 
nitude of the Tfe anomaly decreases again and vanishes 
around 12 T. This disappearance is a consequence of the 
strong hysteresis of the FE transition and occurs only in 
the FC measurements of a&. The upper inset of Fig. ^3 
compares the FC and zero-field-cooled (ZFC) measure- 
ments of ctb for H = 12 T || b. The ZFC curve displays a 
large anomaly at Tfe, which is absent in the FC curve. 
This difference arises from a decreasing T^e(H) with in- 
creasing field. When the lower Tfe (if) is smaller than 
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FIG. 4: (Color online) Field-cooled measurements of ctb(T) in 
fields up to 1 T (a) and between 8 T and 14 T (b). The inset of 
panel (a) compares ctb(T) for H — 4 T applied along b (solid 
line) and c (dashed). The upper inset of panel (b) shows ab(T) 
for H = 12 T (H\\b; dT/dt > 0) obtained in a FC (symbols) 
and a ZFC run (line). The lower inset displays the hysteresis 
between a b (T) for H = 6 T recorded with dT/dt > (solid 
line) and dT/dt < (dashed). Panels (c) and (d) show the 
ZFC magnetostriction obtained for dH/dt > 0. 



our lowest measurement temperature of ~ 4.5 K the FE 
phase cannot be reached in a FC run and, consequently, 
there is no anomaly in the subsequent ab(T) measured 
with increasing T. The FE phase can, however, be en- 
tered by cooling the sample in a lower field, e.g. 6 T (see 
lower inset of Fig. 2t>), or in a ZFC run when the field 
is applied at the lowest temperature. Due to the large 
hysteresis the sample remains in the FE phase up to high 
fields and therefore the ZFC a;,(T) curves show anoma- 
lies at the upper Tfe (if) when the FE phase is left upon 
heating. 

The lower panels (c) and (d) of Fig.0|show the magne- 
tostriction. In agreement with the negative signs of the 
ctb anomalies we observe a pronounced contraction at the 
field-induced ICAFM-to-cAFM transition for T > 13 K. 
At lower T we still find clear anomalies, however, of op- 
posite sign. This sign change results from the fact that 
below about 10 K the FE phase is entered, which has a 
significantly longer b axis than both, the cAFM phase and 
a low-T extrapolation of the ICAFM phase (see Fig- EJ - 
The strong decrease of Ab/b at 10 K for H > 10 T signals 
a field-induced FE-to-cAFM transition due to the above- 
mentioned negative slope of the Tfe (if) boundary. 

Based on our measurements of on and ALi(H)/ Li with 
H\\i for all three crystallographic axes i we derive the 
phase diagrams presented in Fig. |SJ The weak field de- 
pendence of Tn — 42 K well agrees with previous results A 
The new feature of our phase diagram is the down- 
bending of the ICAFM-to-cAFM phase boundary at low 
magnetic field. This transition exhibits a strong hys- 
teresis and we conclude that the ICAFM and the cAFM 
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FIG. 5: (Color online) Phase diagrams of GdMn03 based on 
the measurements of a, (triangles) and ALi (H) /Li (squares) 
measured along i — a, b, and c with H\\i. Filled and open 
symbols are obtained with increasing and decreasing T or H, 
respectively. Hatched areas signal regions of strong hysteresis, 
where different phases coexist. 



phases coexist in the low-field low-temperature region. 
Such a coexistence can naturally explain why the zero- 
field ati(T) curves only show some broad anomalous fea- 
tures around 18 K instead of the distinct anomalies which 
signal the transition from the cAFM to the ICAFM phase 
for larger fields. This conclusion is also supported by 
measurements of the polarization showing that the mag- 
nitude of P in this hysteretic region is much smaller than 
both, P for larger H\\b and P of i?Mn0 3 with R = Tb 
and Dyii 

As shown in Fig. 0Ji the magnitudes of the a& anoma- 
lies at Tfe and T c simultaneously evolve between and 
1 T. This correlation suggests that the ICAFM-to-cAFM 



transition is a precondition for the FE ordering. As 
mentioned above the notation 'cAFM' should be treated 
with some caution, because the magnetic structure of 
GdMnC>3 has not yet been unambiguously determined. 
In a simplified picture, the proposed structure can be de- 
scribed as follows: the Mn moments are oriented approx- 
imately along b with some canting towards c, along the a 
(c) direction neighboring moments are essentially parallel 
(antiparallel) with respect to each other, and along b an 
incommensurate modulation of the moments is present 
above and vanishes below T c . This view is supported by 
the phase diagram, since the stabilization of the cAFM 
phase is most pronounced for H\\c as it is expected for a 
cAFM phase with a weak ferromagnetic moment pointing 
along c already in zero field. A field along a also points 
approximately perpendicular to the Mn moments, and 
one can therefore expect that the Mn moments (and also 
the weak ferromagnetic moment) are slightly canted to- 
wards a. However, a more drastic change is expected for 
H\\b, since 50% of the Mn moments are oriented roughly 
antiparallel to H. In simple antiferromagnets this config- 
uration usually leads to a spin-flop transition, i.e. the Mn 
moment would jump to the ac plane and cant towards b. 
Such a configuration would be analogous to those for H\\a 
or c and according to Ref . IeH none of these configurations 
would lead to a finite polarization via a coupling of the 
magnetic and ferroelectric order parameters. However, in 
view of the more complex magnetic structures observed 
in neighboring i?MnC>3, one may speculate that a simi- 
lar complex structure could be induced in GdMnC>3 for 
H\\b, and this might explain the finite FE polarization via 
a coupling between AFM and FE order parametersiiSii^ 
Hence, there is need for a detailed determination of the 
magnetic structure of GdMnC>3. 

In summary, we have presented a study of the 
magnetic-field temperature phase diagram of GdMn0 3 
via thermal expansion and magnetostriction measure- 
ments. We find that both, the ICAFM-to-cAFM as well 
as the FE transition are of first-order type and strongly 
hysteretic. The hysteresis is most pronounced in the low- 
field range. We find a down-bending of the ICAFM- 
to-cAFM phase boundary and evidence for coexisting 
ICAFM and cAFM phases in this low-field range. 
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